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Context. Several spiral galaxies, as beautifully exhibited by the case of NGC 6946, display a prominent large-scale spiral structure 
in their gaseous outer disk. Such structure is often thought to pose a dynamical puzzle, because grand-design spiral structure is 
traditionally interpreted as the result of density waves carried mostly in the stellar disk. 

Aims. Here we argue that the outer spiral arms in the cold gas outside the bright optical disk actually have a natural interpretation 
as the manifestation of the mechanism that excites grand-design spiral structure in the main, star-dominated body of the disk: the 
excitation is driven by angular momentum transport to the outer regions, through trailing density waves outside the corotation circle 
that can penetrate beyond the Outer Lindblad Resonance in the gaseous component of the disk. 

Methods. Because of conservation of the density wave action, these outgoing waves are likely to become more prominent in the outer 
disk and eventually reach non-linear amplitudes. To calculate the desired amplitude profiles, we make use of the theory of dispersive 
waves. 

Results. If the conditions beyond the optical radius allow for an approximate treatment in terms of a linear theory, we show that fitting 
the observed amplitude profiles leads to a quantitative test on the density of the disk material and thus on the dark matter distribution 
in the outer parts of the galaxy. 

Conclusions. This study is thus of interest to the general problem of the disk-halo decomposition of rotation curves. 
Key words, galaxies: spiral - galaxies: structure - galaxies: halos - galaxies: kinematics and dynamics 
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1. Introduction 



Deep HI observations of nearby galaxies have led to the dis- 
covery of a number of important phenomena that are changing 
our views on the structure and dynamics of galaxies. These in- 
clude the existence in early-type galaxies of regular and radi- 
ally extended HI disks (Oosterloo et al. 2007a), the presence in 
spiral galaxies of extraplanar gas characterized by slow rotation 
(for NGC 891 see Oosterloo et al. 2007b; for NGC 2403, see 
Fraternali et al. 2002), and the properties of small-scale struc- 
tures in the HI distribution (for NGC 6946, see Boomsma et al. 
2008). One interesting related discovery is the existence of reg- 
ular and prominent spiral arms in the gaseous outer disk, well 
outside the bright optical disk (Shostak & van der Kruit 1984; 
Dickey et al. 1990; Kamphuis 1993; for NGC 2915, see Meurer 
et al. 1996; for NGC 3741, see Begum et al. 2005). In this re- 
spect, a particularly impressive example is given by the case of 
NGC 6946 (Boomsma, 2007; Boomsma et al. 2008; see Fig. 1), 
where a spectacular set of gaseous arms can be traced all the 
way out, with a significant degree of regularity and symmetry 
even if the outer disk is clearly lopsided and characterized by a 
fragmentary structure; the outer arms also appear to contain stars 
(see Ferguson et al. 1998, who also analyze the interesting cases 
of NGC 628 and NGC 1058, and Sancisi et al. 2008). 

The study of spiral structure in galaxies has received great at- 
tention in the past. It is now generally thought that grand-design 
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structure is the manifestation of density waves, mostly carried 
by the stellar component of galaxy disks. A general framework 
for the interpretation of the observed morphologies, based on a 
density wave theory, has received significant support from the 
observations of spiral galaxies in the near infrared (see Bertin & 
Lin 1996 and references therein). It might thus appear as a puz- 
zle to find well-organized spiral patterns in the outer disk, in a 
region where stars are practically absent. 

The above-mentioned deep HI studies also serve as interest- 
ing probes in view of defining an appropriate visible matter - 
dark halo decomposition of the gravitational field in galaxies. In 
particular, the studies of prominent spiral arms in the gaseous 
outer disks have raised two issues that separately point to the 
question of whether the outer disk is light or heavy. On the one 
hand, concerns have been raised about the applicability, in such 
region, of the criterion for the onset of star formation proposed 
by Kennicutt (1989), which relies on a threshold on the ax- 
isymmetric stability parameter Q = ck/jtGo-. Accordingly, light 
disks, with low <x, should be unable to make new stars; but in 
the outer parts the disk might be flared and three-dimensional 
effects may change the picture significantly (in this regard, see 
also the comments made by Ferguson et al. 1998). On the other 
hand, there is widespread belief (see Sancisi et al. 2008, p. 212) 
that light disks should be unable to support spiral structure (see 
Toomre 1981; Athanassoula et al. 1987; criticism against this 
belief can be found, e.g., in the article by Bertin et al. 1989a). 

In general, the arguments often put forward in favor of a 
maximum-disk decomposition of the rotation curves of spiral 
galaxies (starting with van Albada & Sancisi 1986) still await a 
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decisive measurement to remove the remaining degeneracy that 
characterizes such a decomposition. Some projects, such as the 
"Disk Mass Project" (Verheijen et al. 2004, 2007), aim at making 
full use of three-dimensional gas and stellar dynamical data so 
as to decompose the field, much like in the classical problem of 
the disk thickness in the solar neighborhood (Oort 1932, Bahcall 
1984, Kuijken & Gilmore 1989, Creze et al. 1998, Holmberg & 
Flynn 2000, 2004). Dark halos are generally thought to be made 
of collisionless dark matter and to have spheroidal shape, but 
it would be important to have direct proof that the outer disk 
of spiral galaxies is indeed light, in contrast with the possibility 
that the disk be heavy because of large amounts of molecular 
material (Pfenniger et al. 1994; see also Revaz et al. 2009). 

One interesting aspect of galactic dynamics is that models 
for the interpretation of observed structures generally offer an 
independent diagnostics of the overall mass distribution. In this 
paper we will present one more example of this general aspect of 
dynamics. Indeed, we will propose a model for the observed spi- 
ral structure in the outer disk and then will show how the model 
can be tested and applied to probe the structure of the outer disk 
in relation with the problem of the amount and distribution of 
dark matter. 

In passing, we note that, in principle, one might test the prop- 
erties of the final model (i.e., disk-dark matter decomposition) 
identified by our technique in specific cases against the expecta- 
tions of non-Newtonian theories of gravity such as MOND; but 
of course, a discussion within MOND of the full problem, in- 
cluding the behavior of density waves, is currently not available. 

The picture presented in this paper is the following. Global 
spiral modes are driven by the transfer of angular momentum 
to the outer regions (see Bertin et al. 1989a,b and references 
therein; see also Lynden-Bell & Kalnajs 1972, Bertin 1983). 
Outside the corotation circle, the transfer is performed by short 
trailing waves. At the Outer Lindblad Resonance, such outgoing 
waves are fully absorbed in the stellar disk (Mark 1971, 1974), 
but only partially absorbed in the gaseous component (Pannatoni 
1983; Haass 1983), so that the signal can penetrate beyond such 
resonance and propagate in the HI disk, if present. The outer spi- 
ral arms are thus interpreted as the natural extension in the outer 
disk of the short trailing waves that are responsible for exciting 
the global spiral structure in the star-dominated main body of 
the galaxy disk. The amplitude profile of such outer arms should 
simply conform to the requirements dictated by the conservation 
of wave action (Shu 1970; Dewar 1972). Because of this conser- 
vation, the amplitude of the outer arms is expected to increase 
with radius, in the regions where the inertia of the disk becomes 
smaller and smaller, much like ocean waves can reach high am- 
plitudes when moving close to the shore. In these outer regions, 
the density wave is thus carried by the gas, but the stars present 
would collectively respond and some new stars may be born be- 
cause of gas compression, following the gaseous arms. 

Of course, galaxies such as NGC 6946 or NGC 628 and the 
Blue Compact Dwarf NGC 2915 are very different objects; each 
case should thus be studied separately in detail and each individ- 
ual object may have its own special character. Here we wish to 
offer one quantitative reference frame for a common mechanism 
that in general should operate in the outer parts of galaxy disks. 

A fully non-linear, three-dimensional analysis of the proper- 
ties of density waves in the HI outer disk is not available, but 
we can hope that under suitable conditions an approximate de- 
scription based on the laws of conservation of wave action for 
low-amplitude density waves is viable. In any case, it should be 
tested against the observations. Since fitting the data requires as- 
sumptions on the density associated with the spiral arms and the 




Fig. 1. Total HI distribution of NGC 6946 superimposed to its 
optical image; courtesy of Filippo Fraternali (see Boomsma et 
al.2008). 

density associated with the fluid basic state, a test of this sce- 
nario would be able to tell how much mass in the outer disk is 
present in the form of molecular gas. 

2. Physical model and calculation of the amplitude 
profiles 

We refer to a grand-design spiral galaxy dominated by a global 
mode with m arms and pattern frequency Q. p and, for simplicity, 
consider a barotropic fluid model of an infinitesimally thin disk 
as an idealized representation of the gaseous outer regions of the 
galaxy. In terms of standard polar cylindrical coordinates (r, 9), 
let cr = o"o + <T\ be the disk mass density, (rQ + v{fi 'g + it\~i ' r the 
fluid velocity field, and c the effective sound speed of the fluid. 
The quantity Q. = £l(r) is the differential rotation. As a measure 
of the distance from the corotation radius, we then refer to the 
dimensionless quantity v = m(Q p -Q)//f, where k is the epicyclic 
frequency. In the following, subscript refers to the axisymmet- 
ric basic state of the disk and subscript 1 to the perturbation, 
on such equilibrium, associated with the large-scale spiral struc- 
ture. In the notation just introduced, the relevant axisymmetric 
stability parameter is defined as Q — ck/(7tGcto)- 

According to the general picture of the modal theory of spiral 
structure (see Bertin & Lin 1996, and references therein), outside 
the corotation circle the global mode is associated with an out- 
going (short) trailing wave. In general, for a normal spiral mode 
the corotation circle is expected to occur at the edge of the opti- 
cal disk (i.e., at r co « 3 - Ah, with h the exponential scale of the 
stellar disk), while for a barred spiral mode the corotation circle 
is expected to occur just outside the tip of the bar, in the middle 
of the optical disk, at r co * 2h. In the gaseous outer disk, the 
outgoing trailing wave can penetrate beyond the Outer Lindblad 
Resonance (which occurs at the radius rout where v = 1). This 
Section is aimed at calculating the amplitude profile of the spiral 
arms in the outer regions, for r > row- In this region the calcu- 
lation of the amplitude profile is particularly simple; instead, the 
properties of spiral structure out to r ~ row are determined by 
the processes that govern the main body of the disk (see Bertin 
& Lin 1996). 
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In the linear theory of density waves, the calculation can be 
performed in a straightforward manner by imposing the conser- 
vation of the density of wave action (see Shu 1970; Dewar 1972; 
Bertin 1983). Mathematically, this is equivalent to carrying out 
the analysis that leads to the standard algebraic dispersion rela- 
tion for density waves D(v, \k\) = 0, with D(v, \k\) — v 2 - 1 - 
Q 2 k 2 /4 + \k\ (where the radial wavenumber k = Ick 2 /{IjtGctq) 
is associated with the pitch angle of the spiral arms, tan/ = 
m/(rk)), to the next order in the WKB expansion (e.g., see Bertin 
2000, Chapter 17.3; see Eqs. (17.2), (17.22), (17.28), (17.30), 
and (17.42)). 

In particular, the linear density wave analysis of a zero- 
thickness barotropic fluid disk is best carried out in terms of the 
enthalpy perturbation hi = (c 2 Ictq)o-\. The WKB asymptotic 
analysis then leads to a Schroedinger-like equation of the form 



we find: 



u" + g(r, u))u — . 



(1) 



where, for tightly wound density waves, the function g is given 
approximately by 



g(r, oS) = 



- 1 + V 



(2) 



The function u in Eq. (Q]) is defined in terms of the function h\ in 
such a way that the amplitudes of the two functions are related 

as 



Ifcil 2 = 



rcr a 



(3) 



The WKB analysis of the turning-point equation for u shows 
that, away from the turning points, the amplitude of u must scale 
as \u\ ~ Igr 1 ^ 4 , which is interpreted as the conservation of the 
density of wave action. A little algebra then shows that such re- 
lation is equivalent to the condition 



1± 
where 



cc G(v, GKWo 4 , 



G(v,Q) = 



v l -\ 



e 2 Vi + (v 2 -i)e 2 



(4) 



(5) 



The proportionality constant is independent of r. Thus Eq. dU 
provides the desired expression for the amplitude profile of the 
density wave in the outer regions. 

In order to complete the description of the amplitude profiles 
of the spiral arms, we may then consider the linearized continu- 
ity equation: 



3o "l 1 9 , n 1 9 , c\\ t\ 

+ "7-(r<7-0«i) + -"^(covi + o-irLl) = . 
at r or r oO 



(6) 



Since the vorticity equation shows that vi = (i / v)[k / (2Q.)]u\, so 
that i'i and u\ have the same order of magnitude, in the continuity 
equation we can neglect the vi contribution consistent with the 
present WKB approximation, thus obtaining: 



(TokU\ ~ VKCTl . 



(7) 



By inserting here the solution for k associated with the (short) 
trailing wave branch, 



(8) 





l K 




0"1 




— v — 
\Q 


)[\rk\) 


O"0 



(£K' Q)r 



-3/2-1 
V n 



with 



H(v,Q) = 



vQ 2 VG(^0 

i + vi + (v 2 - i)e 2 



(9) 



(10) 



Again, the proportionality factor is independent of r. Note that 
the first part of Eq. (0 is an equality, not a proportionality re- 
lation. In other words, in the context of the linear WKB theory, 
the scale in velocity amplitude u\ is uniquely determined by the 
scale in the amplitude of the density wave cr\. Note also that 
the perturbation Mi, together with v\, will generate the "wiggles" 
in the velocity field that are characteristic of density waves (for 
M81, see Visser 1977); in the observed cases where the ampli- 
tudes are large, a non-linear theory is required for a quantitative 
comparison with the observations. 

2.1. Effects of the flaring of the outer disk 

The analysis described previously can be generalized to include 
the effects of finite thickness of the disk. In the discussion of 
the dynamics of self-gravitating disks, these effects are often ig- 
nored, but may actually play an important role; in our case, sig- 
nificant effects would be naturally expected if in the outer parts 
the disk is flared, as often observed. Qualitatively, these effects 
should become significant for short waves, when the wavelength 
of the density wave becomes comparable to the thickness of the 
disk; in terms of local stability, finite-thickness effects are known 
to be stabilizing, because they effectively "dilute" the gravity 
field. Quantitatively, they are approximately incorporated by the 
following dispersion relation (see Vandervoort 1970, Yue 1982): 



D ft {v, \hzo) = v 2 - 1 - l -Q 2 k 2 + -^X- 
4 1 + \k\zo 



= 



(11) 



where zo represents the disk thickness (defined in such a way 
that the disk surface density <t is related to the volume density 
po on the equatorial plane by the expression <xo = 2poZo) and 



zo = 



K 2 Zg 



(12) 



The modification of the dispersion relation with respect to 
the standard one, used earlier in Sect. 2, changes the expres- 
sion of the trailing wave-branch that carries angular momentum 
outwards k = k$r = ksTK 2 /(2nG<To), with respect to Eq. ((SJ. 
But the calculation of this wavebranch is straightforward, be- 
cause the new dispersion relation (fTTT i is just a cubic in \k\. Here 

k ST = ksr(v, Q,ld)- 

From the general theory of dispersive waves, we know that 
the flux of density wave action is T — c g Jl, with the group 
velocity c g = -(doj/dk) and the wave action density J?l oc 
cr 2 (dD/d(L>). Thus we have T oc (dD/dk). In the derivation re- 
ported in the first part of the Section the factor (dD/dk) enters in 
the conservation equation ©, as (dD/dk) = yi + (v 2 - l)Q 2 , 
with the latter derivative evaluated from the zero-thickness dis- 
persion relation on the short trailing wave-branch ([8]). 

Therefore, the desired conservation equation in the finite 
thickness case becomes: 



G f ,(v,Q,zo)r l K 4 cr Q 4 , 



(13) 
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where now we have 



Gf t (v, Q, zo) = 



v l -\ 



Q 2 (dD fl /dk) 



(14) 



here the partial derivative is evaluated at k — ksr- 

The corresponding expression for the amplitude profile in the 
velocity field is obtained from Eq. (0 by replacing H(v, Q) with 
Hf t (y, Q,zo), with 



H f ,(v,Q,zo) = 



V-y/G/rO, Q,Zo) 

fer(v,e,z )l 



(15) 



In the following subsection it will be shown that flaring ef- 
fects are under control and do not change significantly the gen- 
eral predictions of the zero-thickness theory. 

2.2. A simple reference model 

To test the overall picture we may consider the following simple 
reference model. We refer to a two-armed spiral structure (m = 
2) in an outer disk characterized by a flat rotation curve, so that 
k = V2f2, Q/Op = r co /r, and v = V2(r/r co - 1). To be specific, 
we will assume that r co = 3/z, in terms of the exponential scale 
h of the (inner) stellar disk. Thus the Outer Lindblad Resonance 
will occur at rout ~ 5.\2h. Therefore, at r,-„ = 6h we are outside 
the circle associated with the Outer Lindblad Resonance, in a 
region where we expect the disk to be gas dominated, so that 
beyond such radius we may proceed to apply our fluid model. 

As a further simplification, we take the conservative case in 
which cr = <x (r,„)(6/!/r). This corresponds to a very gentle 
decline of the gas density profile. (For faster declining profiles, 
we expect that the relative strength of the spiral amplitude should 
be more pronounced as we move outwards (cf. Eq. ©), and thus 
give rise to a stronger effect.) In our simple reference model, the 
gas density is then proportional to the differential rotation, so 
that we find |cri/cr | 2 oc G(v,0/r and |«i/(rQ)| cc H(v, Q)/r l/2 . 

In a zero-thickness disk, within the adopted reference model, 
the pitch angle of the spiral arms depends only on v and Q, be- 
cause (cf. Eq. ©) 



| rfe | K i_(l+ Vl + (r 2 - 1)Q 2 ) , 



(16) 



where the proportionality constant is independent of radius. In 
this case we are thus left to discuss the dependence Q = Q(r). 
There are observational indications (cf. Boomsma et al. 2008, 
Fig. 6) that the gas velocity dispersion is steadily declining in the 
outer disk, while in the present simple model the quantity k/ctq 
is r-independent. So we might be led to argue that a declining 
g-profile would be realistic. However, the observed decline in 
the gas velocity dispersion is just likely to reflect the fact that, 
because of the stabilizing role of thickness (see previous subsec- 
tion), the disk can get colder and colder (in terms of the standard 
Q) and still remain on the margin of local axisymmetric insta- 
bility. Then, for the present zero-thickness reference model we 
think it appropriate to consider Q « constant, and, for simplic- 
ity, we take Q — 1 . Note that this condition of marginal stability 
is not strictly necessary, because in the outer disk short trailing 
waves can propagate even at higher values of Q. 

In conclusion, the present simple zero-thickness model is 

characterized by G(v,0 = (v 2 -l)/v, H(v, Q) = ^v(v 2 - 1)/(1 + 
v), and rk oc (1+v). By applying Eqs. (0, (|4j, and (O we can pro- 
ceed to calculate the desired profiles z'(r), |o"i/o"o|, and |«i/(rQ)|. 




Fig. 2. For our simple reference fluid model, with ctq oc 1 /r and 
Q oc 1/r, the thick solid line represents the marginal stability 
condition Q = Q max {zo) and the dashed curve the corresponding 
condition for the case in which the dilution of the gravity field 
(associated with the finite thickness of the disk) in the relevant 
dispersion relation is described by an exponential factor (instead 
of the rational factor used in Eq. (fTTI) ). The thin rising curves 
represent the function Q = QfUudizo) for the two cases of a fully 
self-gravitating layer (upper curve) and for a non-self-gravitating 
layer (lower curve). 



In the finite thickness case, the wavenumber for short trail- 
ing waves ksr(v, Q,z.o) is obtained from the dispersion relation 
Eq. (II It . which can also be used to calculate the two relevant 
quantities G/,(v, Q,zo) and Hf,(v, Q,zo). The remaining point 
that requires discussion is the radial dependence of the two func- 
tions zo = zo(r) and Q = Q(r). 

As to the vertical thickness of the outer disk, for our refer- 
ence model, characterized by Q. ~ 1/r and cr ~ 1/r, it can be 
shown (see Appendix A in Bertin & Lodato 1999) that, in each 
of the two opposite limits of a totally self-gravitating disk and 
of a non-self-gravitating layer, the thickness behaves as zo ~ >*, 
so that taking £0 ^constant is a reasonable choice. In addition, 
since we are referring to a fluid, the relevant velocity dispersion 
tensor is isotropic (c r — c z — c), so that there is a one-to-one 
relation between the value of Q and the value of zo that we are 
going to take; let us denote this relation by Q = Qfi„id(zo). On 
the other hand, for a proper comparison with the zero-thickness 
reference model described earlier in this Subsection, we should 
assume that the disk is at marginal stability, i.e. Q = Q ma x(zo) 
(see Eq. (15.22) in Bertin 2000). By combining the above re- 
quirements into QfMdizo) = Q = QmaAzo) we get a unique value 
for the pair {zo, Q), as demonstrated in Fig. 2. With this determi- 
nation, from Eqs. (fTTT i. ( TT4l >. (Q3), and (fT3l > we can calculate the 
desired profiles i[f), |cri/o"o|, and |«i/(rQ)| for the finite thick- 
ness model. 

The results of this analysis are illustrated in Fig. 3 in the 
radial interval 6h < r < lOh. In this figure, we have assumed 
that at r = 6I1 the relative amplitude of the density wave is 
o"i /o"o = 0.15 and that the pitch angle of spiral structure at such 
inner location is 15 degrees. The figure shows that even for the 
presently assumed very gentle decline of the gas density distribu- 
tion the amplitude of the spiral wave steadily increases with ra- 
dius. Note that in the two limits, of a totally self-gravitating disk 
and of a non-self-gravitating disk, the finite-thickness effects do 
not change the general picture; similar results and a similar con- 
clusion have been checked to hold using a dispersion relation 
alternative to Eq. (fTTb . in which the dilution of the gravitational 
term tin is exponential rather than rational. 
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Fig. 3. Relative amplitude profiles of the spiral arms in terms 
of density o-i/ctq and of radial velocity u\/(rQ) (dimensionless, 
left axis) and pitch angle i of spiral structure (degrees, right 
axis), for three cases of the simple reference model described 
in Sect. 2.2: zero-thickness disk (thick lines), finite-thickness 
totally self-gravitating disk (dashed lines), finite-thickness non- 
self-gravitating disk (dash-dotted lines). 

3. Discussion and conclusions 

HI observations of the gaseous outer regions provide measure- 
ments of the HI gas density <x H/ , of the rotation curve D., and 
of the turbulent velocity cuu which we may identify with c. We 
may argue that the actual disk density be traced by the atomic 
hydrogen, so that <x = /cr H/ . In the simplest model we may take 
/ ~ 1 .4, i.e. take that the outer disk just contains the relevant pro- 
portion of helium. For the following discussion, we may assume 
that the thickness zo of the gaseous outer disk is not well con- 
strained by the observations, since we wish to consider objects 
with prominent observed spiral arms, and therefore galaxies that 
are not edge-on. On the other hand, instead of considering the 
thickness profile zo = Zo(r) as an additional free function of the 
problem, we may refer to a self-consistent estimate of such pro- 
file, which is readily available (e.g., see Appendix A in Bertin & 
Lodato 1999 and Sect. 2.2 in the present paper). 

Clearly, if the linear theory happens to be viable, the basic 
relations are over-constrained by the data, since we only have 
some leverage on the precise value of Q. p (which then sets the 
form of the function v(r)) and basically no other free parameters. 
Note that the morphology, or a suitable Fourier decomposition of 
the observed structure, would identify the dominant value of m 
and the pitch angle z'(r), i.e., the function k(r). The linear the- 
ory makes specific predictions, such as Eqs. ®, ©, © (or the 
corresponding equations recorded in Sect. 2.1 for the study that 
includes the effects of finite thickness). Much as for the classical 
tests of the density wave theory (for M81, see Visser 1977), we 
may hope that the various observed quantities all fall within a 
reasonable agreement with the theoretical expectations. 

Suppose that we start from such a straightforward compar- 
ison with the observations on the basis of the linear theory de- 
scribed in the first part of Sect. 2. If this attempt turned out to be 
unsatisfactory, we would have three levels of action in interpret- 
ing the data. 

(i) One possibility would be to make use of an HI frac- 
tion / constant, but significantly larger than unity. Physically, 
a choice of this kind would correspond to imagining a heavier 
outer gaseous disk, but still in proportion to the observed HI 
disk, and thus with little relevance to the overall problem of dark 
matter. For given values of cr HI and cuu this would allow us to 



reduce the value of the wavenumber scale k 2 /{IttGctq) and of 
the axisymmetric stability parameter Q, while leaving the rela- 
tive density amplitude /<x unchanged on the left-hand-side of 
Eq. @. Except for a small adjustment through the functions G 
and H, such a constant / would have little or no effect on the fits 
to the observed amplitude profiles dictated by Eqs. (01 and (O. 

(ii) A second possibility would be to make use of a free 
function f = f(r), with the general requirement that f > 1 .4. 
Resorting to such a function would correspond to imagining that 
significant amounts of dark matter would be in the thin disk, pos- 
sibly in molecular form, with a distribution different from that of 
the cold atomic hydrogen. This may thus go in the direction of 
an alternative picture (with respect to the standard picture of a 
spheroidal halo) for the general problem of dark matter (e.g., see 
Pfenniger et al. 1994). Clearly a non-constant HI fraction f(r) 
would change the character of the observed amplitude profiles 
predicted by Eqs. (|4]l and ([§). 

(iii) Finally, it may well be that the model developed above, 
which relies on the predictions of a linear theory, turns out to be 
inadequate. In other words, we should develop a model in which 
the role of the non-linearities associated with the finite ampli- 
tude observed in the prominent spiral arms is properly taken into 
account. 

While a more advanced model, of the kind outlined in item 
(iii) above, is being developed (in this respect, we are encour- 
aged by the fact that a model of the outer disk as a fluid disk is 
likely to be appropriate and that the physical picture considered, 
of an outgoing wavetrain, is relatively simple), in this paper we 
argue that a first test of the simple picture presented here is worth 
trying. If a satisfactory agreement with the linear theory could be 
obtained with fx 1 .4, we would have one additional convincing 
argument that the picture of a spheroidal dark halo indeed holds 
also for the outermost regions of spiral galaxies. In any case, the 
considerable effort required by setting up such a test on a specific 
case, e.g. for NGC 6946, would prepare the ground for a test of 
the more realistic non-linear theory that we plan to investigate 
soon. 

Finally, given the picture provided by the simple reference 
model described in Sect. 2.2, we may argue that for those galax- 
ies in which the outer gaseous disk density declines too sharply, 
the non-linear effects that would rapidly take place because of 
the sharp increase in relative amplitude of the density waves may 
actually tend to break the grand design outer structure and result 
in turbulent dissipation, as often argued in the physical discus- 
sion of the outer boundary condition for the establishment of 
global spiral modes (e.g., see Bertin & Lin 1996, p. 222). 
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